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Abstract. 

Faraday's and Furry's formulae for the electromagnetic momentum of static charge 
distributions combined with steady electric current distributions are generalised in order to 
obtain full agreement with Poynting's formula in the case where all fields are of class C^, i.e., 
continuous and continuously differentiable, and the integration volume is of finite or infinite 
extent. 

These three formulae are further generalised to the case where singularities are allowed 
to exist at isolated points in the fields, and at surfaces separating domains in which the 
distributions are of class C^. 

Applications are made to electric and magnetic, point-like and finite dipolar systems, 
with an emphasis on the impact of singularities on the magnitude and location of the 
electromagnetic momentum. 

Published: Eur. J. Phys. 28 (2007) 1021-1042. 
1. Introduction 

As recently recalled by Jackson [[U and McDonald [2] , for systems in which effects of radiation 
and of retardation can be ignored, the total electromagnetic momentum can be calculated in 
various equivalent ways, 

^Poynting = JJJ ^^V^^, (1.1) 
Ppurry = jjj ^V^, (1.2) 
^^Faraday = jjj d^V^, (1.3) 

where p is the electric charge density, A is the magnetic vector potential (in the Coulomb 
gauge where V ■ A = 0), E is the electric field, B is the magnetic field strength, (p is the 
electric (scalar) potential, j is the electric current density, and c is the speed of light (which 
we take as c = 1 in the following). Finally, it is generally assumed that these formulae apply 
when the integrals extend over the whole three-dimensional space, and when all the integrands 
decrease sufficiently fast at infinity. 
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In this paper we show that the electromagnetic momentum of a finite or point dipole 
in an electromagnetic field can be consistently calculated with the three equivalent formulae 
(ll.H - ll.3l) . as well with their more general forms applicable when the integration volume is 
finite, using the methods presented in reference [3Jg| These methods allow a straightforward 
calculation in all cases, including those in which there are singularities in the fields. 

The interest of calculating the electromagnetic momentum of a stationary dipole in 
a constant electromagnetic field is not just to confirm that such systems have a "hidden" 
electromagnetic momentum even though they are static charge distributions combined with 
steady electric currents [|4l[5]|: these investigations confirm that the electromagnetic fields of 
point singularities have to be supplemented by (5-like contributions which cannot be neglected 
in the general case E |7l [H |3 • 

For example, the complete field strength of a point magnetic dipole of magnetic moment 
m is ll3l 

^^(-1 = (3^7r^ - ^)^^^) + — -2 — (1-4) 

where u is the unit vector in the direction of r, and, as will be recalled in section [3l the 
generalised function T insures that differentiation at the position of the singularities properly 
leads to the 5 functions which arise when calculating the fields and currents according to 
Maxwell's equations. For instance, after volume integration, the 5 singularity in (11.41) gives 
the (5-like term discussed by Jackson in IflOl p. 184], which is essential in calculating the 
hyperfine splitting of atomic states [|7|. 

Similarly, the complete field of a point electric dipole (which is also called electric 
'dimonopole' to stress that its distributional properties are radically different from those of an 
intrinsic magnetic 'dipole') of moment p = qd, where d is the infinitesimal vector separating 
the positions of two poles of charge ±q, is O 

EUr) = [^^, ^)^(^) ^^(0- (1-5) 

Here we have the 5 singularity that after volume integration yields the 5-like term discussed 
by Jackson in IfTOl p. 141], which together with the corresponding one in (11.41 ) explicitly 
appears in the formulation of quantum electrodynamics of point particles having anomalous 
electromagnetic moments HSj. In section [5] it will be seen that this 5 singularity is also 
indispensable to calculate the electromagnetic momentum using (11.11 - 11.31) for a point dipole. 
The plan of this paper is as follows: 

In section [2] we derive the Furry and Faraday formulae and keep them in the general 
form including the surface terms. It is assumed that all potentials, fields, and source charge 
or current distributions are well behaved functions, i.e., differentiable with a continuously 
differentiable derivative at each point. 

In section [3] we consider the case where the potentials, fields, and sources may have 
singularities, such as those at the locations of point charges or dipoles, as well as those at the 
boundaries of idealised finite systems. 

I There are a few misprints in that paper which are corrected in the errata given at the end of the present paper. 
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In section |4] we calculate the electromagnetic momentum of a point electric dipole in a 
locally uniform magnetic field Bq. It will be found that the dipole's momentum is Bq x p/3, 
provided the surface terms are included in the Faraday and Furry formulae, and provided the 
5-function term appearing in the full expression of the electric dipolar field, i.e., equation (|1.5I) . 
is included when evaluating Poynting's expression. 

In section!!] we calculate the electromagnetic momentum of a point magnetic dipole in the 
field of a point electric charge. It will be found with all three methods that the total momentum 
is Eq X m, where Eq is the electric field at the position of the dipole, provided the (5-function 
term appearing in equation (11.41) is included when evaluating Poynting's expression. 

In section |6] a finite electric and magnetic dipolar configuration is considered: two 
concentric shells, one "dipole charged" and the other "dipole magnetised" [fTTTl . Since such 
shells produce a field that is dipolar outside them and uniform inside, this configuration 
provides an idealised model for a finite electric dipole locate within a uniform magnetic field, 
or vice versa, which has the advantage to be exactly solvable [5J. It will be found with Furry 's 
and Faraday's methods that the total electromagnetic momentum is equal to Bq x p/2. This 
result will be confirmed with Poynting's method, which will also enable us to find out where 
the electromagnetic momentum is located, as well as the respective amounts of the various 
regular (i.e., field-like) and singular (i.e., particle-like) contributions to this momentum. 

Finally, in section |71 we come back on our main results and end with a brief discussion 
of the broader implications of the results obtained by applying Poynting's method to the 
double- shell model of section [61 which show that the integrated electromagnetic momentum 
can be considered as being entirely located in the fields, or entirely bound to the sources, or 
else partially located in the fields and partially bound to the sources. 

2. Generalised Furry and Faraday formulae 

In this section we derive Furry 's and Faraday's formulae in their general form starting from 
Poynting's formula (11.11) . i.e., 

i^Poynting = JJJ ^'V^^. (2.1) 

This assumes that Poynting's formula is the fundamental expression of the electromagnetic 
momentum, which is generally accepted to be the case in the relativistic theory, and which 
is taken as the starting point in Jackson's derivation of these formulae [IJ. However, as can 
be seen for example in Furry's paper [4J, there is much flexibility to move between the three 
formulae, provided some supplementary assumptions are made. It is in order to clarify these 
assumptions that we rederive Furry's and Faraday's formulae here. 

Also, we would like to clarify to what extent the existence of singularities in the integration 
volume over which Poynting's formula is evaluated affects the validity of Furry 's and Faraday's 
formulae. In view of this we will first derive them assuming that there are no singularities, 
and consider the case where there are singularities in section [3l 

As will be seen, the main tool for deriving Furry's and Faraday's formulae from Poynting's 
(and vice versa) is Gauss's theorem. In order to make integrations by parts and use that theorem 
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to transforai volume integrals into surface integrals it is necessary that the integrands should 
be sufficiently smooth functions, i.e., of class C^, that is differentiable with a continuously 
differentiable derivative at each point of the volume. We will therefore agree that the symbols 
V and d^V correspond to a simply connected volume, and that the integrands are finite and of 
class at all points within it. 

Concerning the validity of all three formulae there is also the question of the convergence 
of the integrals, especially when the integration volume is infinite. As the answer to this 
depends on the specific integrand, we will for simplicity give only a sufficient condition, based 
on the demand that all functions are "decreasing sufficiently fast at infinity." For instance, the 
criterion 

1^1 |S|(r ^ oo) = C(r-^-^), (2.2) 

where < e < 1, insures that the integral in Poynting's formula behaves as 0{r^'') at infinity. 
Similar criteria will be given for Furry 's and Faraday's formulae to specify sufficient conditions 
under which the surface terms can be neglected. 

2.1. General Furry formula 

Let us express E in terms of the electric potential diS E = — V0. This enables us to rewrite 
dO as 

Ppoynting = jjj >< ^' (2.3) 

where the arrow indicates that the gradient operator acts towards the left. To be able to make 
an integration by parts and use Gauss's theorem we need the gradient to operate on the whole 



Ppurry = jj S xB + HI d'V <Pj . (2.5) 



integrand, that is on both sides of V in the last equation. This is done by rewriting (|2.3I) as 

where the gradients operate as indicated. We can now use Gauss's theorem to replace the first 
volume integral by a surface integral, and Maxwell's equation V x i? = Airj, to get Furry's 
formula in its general form, i.e., 

1 

This expression is equal to Furry's simplified formula (II. 2|) when the surface term can be 
neglected, that is when 

10 B\{r ^ oo) = 0{r-^~'), (2.6) 

because the surface element is C(r^). 

2.2. General Faraday formula 

This time we express B in terms of the vector potential as B = V x A, and we use the identity 

ExC^ xA) = "^{A • - ■ '^)A, (2.7) 



On the electromagnetic momentum of static charge and steady current distributions 5 

where the underline indicates the range of the gradient operator. This enable us to rewrite 
(O) as 



fl! SV^{A- E)-^ III £V{E ■ ^)A. 



^Poynting = ^ JJJ ^ V ■ IZ/j - — /// a K ^22/ ■ V (2.8) 

In the first integral we use E = — V0, and we rewrite the second one in terms of a left-right 
operating gradient so we can later use Gauss's theorem. Thus 



Ppoynting = ~ ^ fff d^V^ (A ■ V 



(2.9) 

AnJJJ "^'^^^-^^^^h III d'V{E.^)A. (2 AO) 

For the same reason we introduce a left-right operating gradient in the first integral, i.e., 

^ /// £V^{A-^(P) + ^ III £V^{A-^(P) (2.11) 



^Poynting - jjj • _^ ... , ^^^^^ 

^'^^^ ■ ^^^^ i III ^'^^^ ■ ^^-^^^ 

If we now postulate the Coulomb gauge, i.e., 

V-1=0, (2.13) 

the integral on the right of (12.1 II) is zero because partial differentiations commute. As for the 
integral on the right of (12.121) we use E - ^ = Airp. It comes 

Ppoynting = " ^ ^ d^Vf{A ■ ^0) (2. 14) 

'I^lll "^'^^^ '^^^^ III ^^'^^^ 

where we just need to use Gauss's theorem twice to replace the first two volume integrals by 
surface integrals to get the general form of Faraday's formula, i.e., 

Pparaday = / [^ '^(^ " S) + {E ■ d^S)!^ + jjj d'VpA. (2.16) 



This expression is equal to Faraday's simplified expression (11.31) when both surface terms can 
be neglected, that is when 

\(j)drAr\{r ^ oo) = 0{r-^''), \E\\A\{r ^ oo) = 0{r~'^-'), (2.17) 

where dr-Ar is the partial derivative of the radial component of A. As for the Coulomb gauge, 
equation (12.131 ). it is a necessary condition for the applicability of both the general (12.161) and 
the simplified (11.31) forms of Faraday's formula. 



3. Inclusion of singularities 

The simplified expressions (|l.lj - ll.3l) and the more general formulae l2.5l and l27T6l are rigorously 
applicable to any configuration of static charges and steady currents provided all potentials, 
field strengths, and source charge or current distributions are of class C^. This means that they 
can be directly applied to any engineered electrical system in which all components have a 
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geometry described by smooth finite functions. They can also be applied to systems in which, 
for the purpose of simplifying the calculations, some components are idealised with Dirac 
delta-functions and Heaviside step-functions, in which case spurious singularities may arise 
as artifacts due to these idealisations, but can in general be discarded on the basis of simple 
plausibility arguments. 

On the other hand if these expressions are applied to systems including elementary objects 
such as point-charges and point-dipoles there will be real singularities which in principle 
invalidate the application of these formulae since there will be points at which the integrands 
will no more be finite. In such cases the general method consists of isolating the singularities 
in small regions enclosing them, and of dealing separately with the regular and singular parts 
of the integrals. However, as is well known, the proper way of dealing with the singular parts 
is still an unsolved problem in electrodynamics, and the best that can be done is to use a 
prescription to "regularise" the integrals and possibly "renormalise" the results if so requiredj§| 
This is what we will do in sections |4] and [51 

There is another class of singularities which arise when an idealised configuration such 
as a charged-shell is considered as an elementary system, that is when it is considered as a 
genuine elementary solution of Maxwell's equation, even though it is not believed that such 
a configuration necessarily corresponds to a real physical object (in the sense that a point- 
charged is believed to correspond to an electron). In this case, as will be seen in the example 
studied in section [6l the integrals can be finite because the sources are no more concentrated 
at a point, but there may still be singularities at the boundaries (e.g., at the surface of a shell) 
which have to be physically interpreted. Also, the integrals have to be evaluated separately in 
each domain defined by the boundaries. For the Furry and Faraday formulae this implies that 
the surface terms have to be carefully evaluated at all boundaries to insure that their respective 
contributions are such that the total given by the outermost boundary is physically meaningful. 

Finally, when considering systems in which singularities overlap, or in which elementary 
objects self-interact, there will appear products of distributions which cannot be dealt with 
using standard distribution theory. While new methods exists to evaluate these products [fT2ll . 
degenerate or self-interacting systems will not be considered in the present paper. 

Since the problems associated with singularities are complex and crucially depending on 
the details of the system under consideration, there is no general method to deal with them. 
Nevertheless, some simplification is achieved by characterising the singularities in such a 
way that their origin is made explicit starting from the potentials of the fields. This has the 
advantage that the singularities do not have to be discovered indirectly on a case-by-case basis 
by looking at the details of all expressions and their derivatives, but arise consistently and 
automatically in accord with the principles of the theory of the distributions. 

A method to do that is explained in reference (Si, which results in expressions such 
as equations 11.41 and 11.51 of the introduction. The basis of this method is to make explicit 
the difference between the linear coordinate variables (r in a polar system; x, y, and z in a 

§ This is obviously not very satisfactory, but at present the only way to justify the assumption that all 
distributions considered by electrical engineers are of class C^, which implies that the singular contributions 
stemming from the charge and magnetic moment of all electrons are systematically ignored. 
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Cartesian system) which are continuous, and the modulus \r\ of the position vector (|r| in a 



polar system; a/x^ + + in a Cartesian system) which is discontinuous at r = 0. In a 
polar coordinate system this is done by writing, for any occurrence of the radius vector f, 

r = rT{r)u{e,<p), (3.1) 

where u is the unit vector in the direction of r. Here, the generalised function T(r) has the 
properties 

{undefined for r < 0, 
forr = 0, and ^T{x) = 5{x), (3.2) 

+ 1 for r > 0, 

as well as 

I drT{r)F{r)= j dr F{r), (3.3) 
which is the counter-part of the defining property of the 5 function in polar coordinates, i.e.|il 
dr 5(r)F(r) = F(0), (3.4) 
where, like in (|3.3I) . F e for the purpose of this paper. 

4. Electric dipole in a locally uniform magnetic field 

In this section we consider an electric dipole located in a magnetic field that is uniform over a 
finite volume. That field may be generated by an external system such as a solenoid, or else 
the volume may be assumed to be sufficiently small that the field is constant over it. For the 
electric dipole we take the full expressions for its potential, field, and charge density, namely 

m 

(pdm{r) = ^p- uT{r), (4.1) 
Edmir^ = (3^^ ^'^(0, (4.2) 



3 p- u 
An 

and for the vector potential, field, and current density of the uniform magnetic field the standard 
expressions 

A = ^BoX f, (4.4) 

B = Bo, (4.5) 
/ = 0. (4.6) 

f We adhere to the modern convention that the symbol S corresponds to the Dime measure, i.e., that its defining 
property is {{S\F)) — F{0) independently of the support of S, which is here the half-space [0, oo] rather than 
[— oo, oo] as in the customary definition of Dirac's delta-function. 
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As the current density is identically zero we can immediately conclude that Furry's formula 
without surface terms, i.e., equation (|1.2I) . would give a zero result, which is why we will have 



to use the general formula (|2.5l) instead. We also remark that A satisfies the Coulomb gauge 
so that Faraday's formula is applicable in principle. 

Finally, for the integration volume we take a sphere of radius R, and since we work in 
spherical coordinates we only need a few elementary formulae, i.e., 

(fv = duj dr r^, and d^S = du uR^ , (4.7) 

duj u{p- u) = -^Vi ^iid jj duj u X (px u) = -^P- (4.8) 

We start with Poynting 's formula, equation (11.11) or (12.11) . and replace E and B by (14.21) 
and (|43l) to obtain 

Ppo.„.n, = III rf3y^^ (4.9) 



1 

An 



dr JJ duj'^^(3u{p-u) -pj X Bo (4.10) 

J dr JJ duj 5{r)u{p ■ u) X Bq. (4.11) 

Ignoring the T factor it is seen that the radial integral in equation (14.101) diverges as 1 /r at r = 0, 
and taking T into account leads toT(r)/r = 0/0 which is undefined. On the other hand, using 
(14.81) for integrating over the angles, this equation would be zero if the result from the radial 
integration was finite. Thus, as the radial and angular integrations cannot be interchanged, the 
result is ambiguous. This ambiguity is clearly related to the 1/r divergent factor, which is one 
instance of the fundamental divergencies which plague classical and quantum electrodynamics. 
As is well known, there is no other cure than ignoring these divergences, for example by 
introducing a "covariant" form-factor, and possibly "renormalising" the result afterwards. In 
the present case, a consistent method is to represent the T and 5 functions by the limits e ^ 
of a pair of functions such as, for example, 

2 r 2 e 

T£(r) := -arctan(-), and 5e{r) := — — — -, (4.12) 

TT e Tx + r^ 

and to keep e very small hut finite so that T^(r)/r = 2/TTe at r = 0. In that case the radial 
integration is finite while angular integration gives zero. 

Equation (14.101) is thus zero, and it remains to calculate equation (14.111) . which using (14.81) 
and (13.41) gives the finite contribution 

^Poynting = ^Bq X p, (4.13) 

which has its origin in the 5 singularity of the field strength (14.21) . 

As everything is finite, we can apply the general Furry formula (|2.5|) . and calculate the 
surface term using (14.11) for (p and (14.51) for B. We get 

^ ' ""'3x3 = -— [[ ^p-uduoR^u X Bo, (4.14) 



4vr JJ ^ Att jj R^' 
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that is, after angular integration, 

3 



^Furry, Eq. EJ = X p, (4.15) 



which agrees with (I4.13I ). 

The first surface term in the general Faraday formula (12.161) contains the expression 
"^{A-d^S) which has to be evaluated assuming that the surface element is constant. Using 



(14. 4|) for A we therefore calculate 

^(1- C) = {C-V)A + Cx (V X A) (4.16) 

= {C-V)^Boxf+CxBo (4.17) 

= {C -11)^30 xu + C X Bq. (4.18) 
Thus, replacing the constant vector C by dtoR^u, 

^(1 ■ d'^S) = duR^^u X Bq. (4.19) 
We now use (14.11) for (p and find 

-ill = -1^11 y. udu^wU X 5o 

= X p. (4.20) 

o 

The calculation of the second surface term in (12.161 ) is straightforward: using (14.21 ) for E and 
(1441) for A we get 

= ^ ■ m)5o X m = -ifio X p. (4.21) 

We finally need the volume term in (12.161) . which is also straightforward to calculate: using 
(1431) for p and (|44l ) for A we obtain 

rfVpl= dr jj du ^^^^'^(^)^^o X f = X p. (4.22) 

It remains to add (l430l) . (|43T]) . and (l432l) . i.e., (1/6 - 1/3 + l/2)5o x to get in total 

Faraday, Eq. ^JB = ^' ^"^'^^^ 

in agreement with (14.131) and (14.151) . 

In conclusion, we have found that by including the surface terms the general Furry and 
Faraday formulae give the same electromagnetic momentum as calculated with Poyntings 
formula, provide the (5-function singularity that is present in the dipolar electric field Edm is 
included, whereas the 1/r singularity in the electromagnetic momentum is discarded. 

As an application of the result obtained in this section, we assume that the volume in 
which the magnetic field is constant is the gap between the plates of a capacitor, which it- 
self is filled by a dielectric material such that its total electric polarisation is p: The total 
electromagnetic momentum in the dielectric material is then Bq x p/3. 
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As another application we transpose our result for an electric dipole in a magnetic field 
to the case of a magnetic dipole in a locally uniform electric field. Indeed, if we consider 
Poynting's formula (14.91) and replace Edm by Eq, and Bq by Bd, we see by comparing (11.41) 
and (11.51) that the only difference is in the angular dependence of the 5-function singularity. 
For a magnetic dipole of moment m we have therefore 

2 - 

^Poynting = g^O X (4.24) 

where the factor 2/3 comes from the angular integration with (14. 8|) of the double cross product 
in(fr4l). 



5. Magnetic dipole in the field of a point charge 

We now consider a magnetic dipole located at r = 0, and a point charge positioned at r = a, 
an example considered by Furry in ||4l . For the integration volume we take the whole space, 
and for the potentials, fields, and current densities the complete expressions including all 
singularities [3J. Thus, for the point charge we have 

0p(r) =g— ^T(|f-a|), (5.1) 
|r — a\ 

E,{r) = qlz±.T{\r- a\) - qlzl_5{\r - a|), (5.2) 
\r — aY" \r — 

pp{r) = q } ^.^ S{\f-a\), (5.3) 
47r|r — 

and for the magnetic dipole 

Mrl = ^T(r), (5.4) 

Bd{f) = 3 ^3^ - - T(r) + ^ ^{r), (5.5) 

3 fh X u 
An 

We start with Poynting's formula, equation (11.11) or (12.11) . and replace E and B by (15.21) 
and (15.51 ). Assuming |a| 7^ we have 5{\f— d\)S{r) = so that the product of the two 5 
terms is zero. Under the same assumption the product of the 5{\f— d\) and T(r) terms is 
also zero. Indeed, after recentering the coordinate system at r = a we get as radial integral 
J dr rS{r) = 0. We are thus left with two terms: 

Ppoynting = fff d'V^^^ (5.7) 



Url =^^^W, (5.6) 



= — f dr If duo q Z t, T(\f- a|) x ^ ( 3u(p- u) - p\ (5.8) 
4vr Jq JJ \r -a\-^ r \ J 

+ — [ dr [[ du q-Z — ^T(\r — d\) x 6(r)(u x (rh x u)) . (5.9) 
4vr Jq JJ \r - a|^ V / 

Contrary to what happened in the previous section's equation (14.101) . the angular integration in 

(15.81) does not give zero because |r — a| ^ \r + d\. Moreover, there are two singular points to 
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consider: f = a and r = 0. Fortunately, the point f = a gives no trouble because the volume 
element in spherical coordinates centered at this point has a factor |r — which cancels the 
factor |r — a|^^ in Ep. 

As for the 1/r divergence, it can be taken care of in the same way as was done in the 
previous section with equation (14.101) . The reason is that very close to r = the difference 
between |r — a| and |r + a| becomes negligible and the net result of angular integration in a 
region close to r = is to yield a zero provided the radial integral over that region is finite. 
Thus, if we replace the zero lower bound in (|5.8I) by a small value e we get a finite result, 
which moreover is independent of the precise value of e. Unfortunately the evaluation of this 
integral is tedious, so that we just quote the result of reference [i4i p. 624-625], i.e., 

1 q 

^Poynting,Eq. HJ = 3^"^ >< (5.10) 

On the other hand, the singular contribution (15.91 ) is easy to calculate because 5{r) implies that 
we can replace r — a by —a. Thus, 

4oynting,Eq. m = [ ^^^^ JJ duj Q-^ X (u X {m X u) 

2 q 

= — -m X a. (5.11) 

3 a'^ 

This partial result has a simple interpretation: It is the electromagnetic momentum of a 
magnetic dipole at a point where the electric field is Eq = —qa/a^, i.e., the electric field at the 
origin in our configuration. Thus equation (15.1 II) agrees with (14.241) . 

Adding (15.101) and (15.111) we obtain the total electromagnetic momentum of a system 
comprising a point magnetic dipole and a point charge, i.e., 
~* Q 

-Ppoynting = X O (5.12) 

This result is easily confirmed by means of Furry 's and Faraday's formulae in the simple forms 
(11.11) and (11.21) because the integrands in the surface terms fall of more rapidly than r^^ for 
r ^ oo, and because the vector potential (|5.4I) satisfies the Coulomb gauge. 
For instance, Furry's formula is 

PP.„. = III ■'"''^^ = III ''"''"'p^Tdf - a\)l'^6(r), (5.13) 

which after developing l/\f—d\ around r = as 
1 1 f ■ a. 



, -(1 + ^), (5.14) 

|r — a| a 

gives a result equal to (15.121) . Similarly, Faraday's formula is 

^ III SV,A ^ III i^V.'^^^^r), (5.,5) 

which after the change of integration variable r — > |r — a| also gives a total electromagnetic 
momentum equal to (15.121) . We have therefore confirmed the total (15.121) by means of Furry's 
and Faraday's formulae, which by subtracting the partial result (|5.1 II) enables us to confirm 
the correctness of Furry's independent calculation of the first partial result (|5.10|) . 
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It is now worthwhile to compare the method used here to that of Furry [4J : In our case 
we have assumed that the dipole was strictly point-like so that we had to use a regularisation 
process typical of quantum electrodynamics to discard the unphysical 1/r singularity in (15.81) . 
while the physical 5 singularity in (15.51) led to the finite (5-function contribution given by (15.91) . 
On the other hand, Furry assumed that close to r = the point-like dipolar field had to be 
replaced by that of a magnetised sphere of very small radius, which gave a contribution equal 
to that of our 5 singularity. 

Finally, it is of interest to rewrite equation (15.121) in the form 

^Poynting = Eq X ITi, (5.16) 

where Eq = —qa/ is the electric field of the point charge at the position of the magnetic 
dipole. Indeed, this expression immediately generalises to the case where any static charge 
distribution is the sole source of Eq, a result obtained by Aharonov et al., i.e., [fT3l Eq. (4)], 
using a variant of the Furry formula BH Eq. (46)]. Thus, (15.161) is the total electromagnetic 
momentum, i.e., the so-called "hidden momentum," of any configuration comprising a point 
magnetic dipole in a constant electric field. 

6. Finite electric and magnetic dipolar system 

In this section we consider a finite electric and magnetic dipolar system made of two concentric 
shells, one of radius a that is "dipole charged," and another of radius b > a that is "dipole 
magnetised" [11 J. The outer shell has a surface current distribution producing a magnetic 
field that is dipolar outside it and uniform inside, while the inner shell has a surface charge 
distribution producing an electric field that is dipolar outside and uniform inside. While this 
configuration is unrealistic in the sense that it would be very difficult to actually build it 
im, it is perfectly consistent with the laws of electrodynamics, and has the advantage that all 
calculations are relatively simple. Moreover, for the purpose of calculating the electromagnetic 
momentum, this system is similar to a spherical capacitor immersed in a uniform magnetic 
field such as considered in 

6.1. Dipolar fields of charge/current carrying shells 

The fields generated by dipolar distributions on shells provide an example of how the 
distributional methods of section [3] can be used in a case where the singularities are spherical 
rather than punctual, because by symmetry they corresponds to a single point of the radial 
variable r. 

For instance, the scalar potential, field strength, and charge density distribution of a 
charged-shell of radius a producing a constant electric field inside it and a dipolar electric field 
outside are 

4>dm{r,r ^ a) = ^^T(r - a), (6.1) 
Edm{r,r ^ a) = [3 — — — jT(r -a) — — d{r - a), (6.2) 
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Pdm{r,r ^ a) = —S{r - a). (6.3) 

As can be seen, these fields go over to equations (14.11 - 14. 31) as a ^ 0. 

Similarly, the vector potential, field, and current density distribution on a conducting- shell 
of radius h producing a constant magnetic field inside it and a dipolar magnetic field outside 
are 

i,(r,r^6) = ^^T(r-6), (6.4) 
B,{f, r^b)= (3^^ - ^)T(r - b) + U{r - b), (6.5) 

— * 3 TKl X 11 

jd{f,r ^b) = ^5(r - 6), (6.6) 

which reduces to the fields of a point magnetic dipole when 6^0, i.e., (15.4145.61) . 

Equations (|6.H - I6.6I) are the expressions for the fields in the regions r ^ a and r ^ 6 where 
they have a dipolar character. As is well-known, the corresponding surface charge and current 
distributions, equations (16.31) and (16.61) . produce fields that are uniform in the regions r ^ a 
and r ^ 6, i.e., inside the shells To obtain the expressions of these fields including their 
singularities we can proceed the same way as we have done for (16.1116. 61) . that is to related 
these singularities to those of a point source by a substitution of the type r — > i? — r where 
R is the shell radius, because in spherical symmetry such a substitution is equivalent to just 
translating a point along the r axis. 

Thus, to find the corresponding expression for the uniform electric field generated by a 
dipolar surface charge distribution we start from the well-known scalar potential of a uniform 
electric field, namely 0(r) = —EQ-urT(r) where we have included the T function taking care 
of the singularity at |r| = 0. Substituting a — r for r in the argument of T gives the potential 
of the electric field, minus the gradient of this potential yields the electric field strength, etc. 
Thus, remembering that T(a — r)' = —5{a — r), 

4>dm{r, r ^ a) = -Eq ■ urT{a - r), (6.7) 

Edni{r,r ^ a) = EoT{a - r) - u{Eq ■ u) r6{a - r), (6.8) 
3 

pdmir, r ^ a) = -— ■u5{a-r). (6.9) 
47r 

Similarly, the vector potential, field strength, and current density distribution of the 
constant magnetic field inside a shell of radius b producing a dipolar magnetic field outside it 
are obtained from the potential ( 14.41 ) rewritten as A(r) = BqX MrT(r)/2 by substituting b — r 
for the argument of the T function, i.e., 

Ad{r, r ^ 6) = ^Bq x urT{b - r), (6.10) 

Bd{r,r ^b) = BoT{b-r) - x {Bq x u) r6{b-r), (6.11) 

jd{f,r <:b) =+^BoXu6{b-r). (6.12) 

f The proof of this basic resuh of potential theory requires however a development in Legendre Polynomials or 
spherical harmonics, as is done, e.g., for of a magnetised sphere in lilOi p. 195]. 
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It remains to relate Eq to p, as well as Bq to rh. This is done by requiring that the 
potentials (l)drnir,r > a) and (l)dm{'i^,r < a), as well as Aii{f,r > h) and ^^(r, r < h), are 



continuous at i? = a, h, respectively. This gives 




We now interpret the fields (16. 11 - 16.121) as those of a finite spherical electric dipole of 
radius a immersed in a uniform magnetic field produced by an ideal spherical dipolar magnet 
of radius b > a. However, before applying Poyntings, Furry's, or Faraday's formulae, we 
need to verify that their conditions of applicability are satisfied. If we exclude the points on 
the respective shells at r = a and r = b, this is clearly the case for all fields (I6.H - I6.12I) since 
they are finite and of class at all other points. This means that we can divide the whole 
space into three concentric domains and calculate separately the electromagnetic momenta by 
all three methods in each of them. If it turns out that the contributions of the surface terms in 
the general Furry and Faraday formulae cancel each other at the boundaries r = a and r = b, 
the contributions from all three formulae will add up to the same total, provided the volume 
integrals do not give problems. To verify that explicitly it would be necessary to calculate all 
these surface terms, what we will not do because from the form of the fields we expect them 
to cancel, just like we do not expect any contribution from the surface terms at infinity since 
all integrands fall off more rapidly than r^^ for r ^ oo. Referring to the derivation of Furry's 
and Faraday's formulae made in reference [IJ, which implicitly assumes that the fields vanish 
sufficiently fast at infinity, we also remark that the boundary terms at finite r should not give 
any problem. It remains therefore to calculate the volume terms. 

6.2. Furry's expression 

With start with Furry's method because it is the easiest to use since 6 > a in our configuration. 
Thus, using (16.11) for and (16.61) for j, we gelQ 



which becomes, after performing the angular integration using (14.81) and expressing m in terms 
of Bq with (16.141) . 



where we replaced the lower integration bound by b because j is zero for r < b. Since 
T(r — a) = 1 for 6 > a, using (13.41) we get, therefore. 



which is the correct value expected for a macroscopic dipole or capacitor ll2ll5l[Tn. 
f One could equally well have used ( |6.12| i for j. 




(6.15) 




(6.16) 



-Ppurry = -Bq X p, 



(6.17) 
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6.3. Faraday's expression 

Here we use (16.31) for p and, since r ^ b, equation (|6.10l) for A, that is the vector potential 
inside the magnetic shell where 5 is a constant. Therefore^ 

dr JJ duj r^^^6{r - a)-Bo x f T(6 - r), (6.18) 

which, after angular integration, becomes 

1 - 

i^Faraday = -5o X rfr 5(r - o) T (6 - r) , (6.19) 

where we replaced the lower integration bound by a because p is zero for r < a, and the upper 
bound by b since we use (16.101) for A. As T (6 — r) = 1 for b > a, using (13.41) we get, therefore, 

^Faraday = ^Bq X p, (6.20) 

which is again the correct value expected for a macroscopic dipole or capacitor. 



6.4. Poynting's expression 

To evaluate Poynting's expression one has to integrate the cross products of the four fields 
(16.21) . (16. 5k (16.81) , and (16.1 II) , which are not very simple due to the T and 5 singularities at a and 
b. Moreover, Poynting's expression has to be evaluated in the three regions defined by the radii 
0, a, b, and oo. Thus, compared to Faraday's and Furry's methods, the full calculation with 
Poynting's method is laborious. For this reason we relegate the details of these calculations 
to the appendix, section [81 and just collect their results in this section. 

However, as they are very simple to calculate, we first give the respective contributions 
to the electromagnetic momentum in the three regions for the case where the singularities are 
completely ignored, that is, when the integrations are made excluding the T and 5 singularities 
at the radii a and b: 

\bo X p, (6.21) 
0, (6.22) 

\b^ X p. (6.23) 

o 

Thus, the non-zero contributions add up to \Bq x p as they should, and the contribution from 
the fields located between the polarised and magnetised shells is zero. Also, referring to the 
calculation made in section |4] for a point-dipole, we see that the subtotal Bq x p/3 for the 
region < r < 6 corresponds to what was obtained in that section after taking care of a 
divergence at r = 0. This divergence does not arise here because the point-dipole is replaced 
by a dipole-polarised shell of finite radius, which implies that the contribution in the region 
< r < a is finite, and that the contribution for a < r < 6 is zero after angular integration. 

f One could equally well have used (|6.91 l for p. 



Ppoy„ting(< a) = jjj d'V^^ 

^Poynting(a < r < 6) = jjj d^V 

Ppoy„.ng(> b)=lll d'V^ 
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Now, recalling what was said in section[3]about singularities, we remark that ignoring the 
T and 5 singularities at the radii a and h corresponds to the "electrical engineering perspective" 
in which they are considered as spurious. On the other hand, anticipating on the interpretation 
that will be put forward in the discussion section, we can also say that discarding these 
singularities corresponds to the "field perspective" in which the electromagnetic momentum 
is calculated by integrating the Poynting vector over the whole space excluding the points at 
which it is singular, that is the points at which the sources reside. 

Therefore, to find the physical picture corresponding to the opposite perspective, i.e., the 
"particle perspective" in which the singularities in the fields are included, one has to calculate 
the momenta in such a way that the respective contributions from the regular and singular parts 
of the fields can be isolated from each other, and this separately in the three regions ^ r ^ a, 
a ^ r ^ b, and b ^ r ^ oo. 

Consequently, to see what happens when the boundaries a and b are included as upper 
or lower integration bounds, we write them as a ± a and b ± P, where a <^ a and (3 <^ b 
are positive infinitesimals, so that the radial part of the integrations made in the appendix are 
defined as follows: 

pa— a fb—f} POO 

/ dr (...), / dr (...), and / dr (...). (6.24) 

Jo J a+a Jb+P 

Thus, as long as the parameters a and (3 are non-zero these integrals give the contributions 
from the fields ignoring the singularities, and when a ox (3 are set to zero they give these 
contributions including the singularities. More precisely, taking the limits a — > and /? — > 
corresponds to ignoring the singularities, i.e., the discontinuities at a and 6, whereas setting 
a = 0, and/or /3 = 0, corresponds to including the contributions from these discontinuities. 

We now look at the results of these integrations, equations (18.111) . (18.221) . and (|8.31l) of 
the appendix, i.e.. 



^ 1 

^Poynting(0 ^r^a-a) = Boxp -— T(a - r) 

O Qj 

^Poyntmg(a + a^r^b-(3) = BoXp^ T(r - a) 

- Boxpl T(6 - r) 
o 





b-f3 

a+a 
b-13 



(6.25) 



^Poynting(& + j3 ^ T ^ (X)) = Bq X p -^T{r - b) 



(6.26) 
(6.27) 
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6 ■ ' b+i3 

For instance, when both a ^ and P ^ the T functions in (16.25146^27] ) are all equal to 
T(x 7^ 0) = 1, which implies in particular that (|6.26l) is zero. Letting a ^ and /? we 
recover the results (I6.21j - [6l23l) . that is the electromagnetic momenta when the singularities 
are ignored. In this case we can also say that the distributional content of (16.25146^271) . which 
quantifies the discontinuous effects of the singularities at a and b, is ignored. 

On the other hand, if either a or /3 is set to zero, equations (16.25146^271) must be interpreted 
in their distributional sense, which means that the T functions at a or 6 take their point- values 
of T(0) = when a or /5 is zero, as stipulated by (13.21) . In particular, whereas (|6.25l) and 
(I6.27|) show an r-dependence that is continuous through the factors r"^^, equation (|6.26l) is 
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strictly discontinuous and thus gives a non-zero contribution only when the integration bounds 
are exactly equal to a and/or h. 

For example, when a = the electromagnetic momentum given by (16.251) for ^ r ^ a 
is zero: the contribution from the singularity at a cancels the contribution of the integral in 
the interval ^ r < a. At the same time the integral (|6.26l) is no more zero but yields a finite 
contribution of BqX p/3 at the bound r = a, that is, exactly the contribution (16.251) gave when 
the singularity at a was ignored: the electromagnetic momentum which was considered as 
being in the field between and a has been "concentrated" at the singularity at a. 

The same analysis can be repeated for the singularity at r = 6, with the conclusion that 
the effect of including or excluding the contribution of that singularity is to concentrate an 
electromagnetic momentum of i?o x p/6 at r = 6, or to leave it at r > b. 

In conclusion, the total electromagnetic momentum is always Bq x p/2, and the effect 
of including/excluding the singularities is to transfer part or all of that momentum to their 
locations, or to leave it in the fields. If we now simplify (I6.25l - [6l27l) by taking the foregoing 
analysis into account, we are led to the following equations: 

^Poynting(0 ^ r ^ a - a) = Bq X p^T{a), (6.28) 

Ppoynting(a + = Bq X p^(^l - T(«)) (6.29) 

Ppoy„ting(& -P) = 5o X p i (l - T(/3)) , (6.30) 

Ppoynting(fe + /? ^ r ^ Oo) = X p^T{l3). (6.31) 

Here a and P are again small numbers which are either zero or non-zero, but since T(r) is 
discontinuous at zero their purpose is just to give results equivalent to those of integrals (|6.25l - 
16.271) when the singularities are included or not. In particular, the total of these equations 
is always Bq x p/2, and the splitting of (16.261) into (16.291) and (16.301) stresses that these 
contributions arise exclusively from the singularities at a or b. 



7. Discussion 



The primary goal of this paper was to better understand under which conditions the methods of 
Poynting, Furry, and Faraday can be used to calculated the electromagnetic momentum of static 
charge and steady current distributions, either taken as whole systems, or as composite systems 
of finite or infinite size. In order to do that Furry's and Faraday's formulae were generalised 
to the case of a finite volume, which led to generalised formulae including surface terms, 
equations (12.51) and (12.161) . which for that reason do not have the simplicity and elegance 
of Furry's and Faraday's original formulae (|1.2I) and (11.31) . Nevertheless, the generalised 
formulae can, in principle, easily be applied to any system, or part of system, in which the 
potentials, fields, and sources are continuous and continuously differentiable functions — that 
is, to realistic systems which are typically considered by "electrical engineers." 

The second main goal of this paper was to calculate the electromagnetic momentum of 
systems in which the potentials, fields, and sources may have singularities and to use for 
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that purpose the method presented in reference [|3]|. As explained in section [3l that method 
consists of making explicit in the potentials the discontinuity which through two successive 
differentiations leads to the singularity corresponding to the sources. 

In section IH the electromagnetic momentum of a point dipole in a locally uniform 
field was calculated with the generalised Furry and Faraday formulae and found to agree 
with Poynting's formula. As this calculation required discarding a diverging term, the same 
formalism was applied in section [5] to a point dipole in the field of a point charge, and shown 
to be consistent and in agreement with the methods of references BU and [fT3l . In both 
applications the ^-function singularity in the field strength of the point-dipoles was found to 
contribute the finite part of the dipoles's contributions to the total momentum after discarding 
the 1/r divergence. 

In the third application, section [6l the location of electromagnetic momentum in a 
composite system made of two concentric shells modelling a spherical electric dipole of 
radius a immersed in the uniform magnetic field produced by a spherical magnet of radius b 
was investigated. 

In sections [6l2 and [6l3, the total electromagnetic momentum was easily found to be 
equal to Bq x p/2 using either Furry's or Faraday's methods, which could be applied in a 
straightforward manner because the finite size of the model insured that there is no singularity 
at the origin (contrary to the point-dipole systems considered in sections |4] and [5]). This total 
momentum was confirmed in section[6l4 using Poynting's method, which also gave additional 
information on the location of the electromagnetic momentum that is not easily provided by 
Furry's and Faraday's methods — unless the complicated surface terms which cancel each 
others at the boundaries are explicitly calculated. 

In section [6l4, the electromagnetic momenta were calculated separately for the three 
regions (0 ^ r ^ a), (a ^ r ^ 6), and (6 ^ r ^ oo), and their significance analysed. The 
results are presented in detail in equations (16. 25146^271) . and in simplified form in equations 
(16. 2 81 - 16. 3 T1) . In both sets two infinitesimal parameters a and /3 enable us to include or exclude 
the contributions from the singularities in the field strengths at the radii a and b: when a = 
or P = the singularities at r = a or r = 6 are included, whereas when a 7^ or /3 7^ they 
are excluded. This further enables us to summarise the results as four symbolic equations, 
(17.11 - 17.41) . which correspond to the four possible ways the partial results (I6.28I - [63T1) can be 
added up by either including or not the contributions from the singularities at a and b: 

• Field perspective (a ^ 0, (3 ^ 0) : 



The electromagnetic momentum is concentrated at the location of the charge and current 
distributions at a and b. 




(7.1) 




(7.2) 
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• "Faraday's" perspective (a = 0, /3 7^ 0) : 

The electromagnetic momentum is at the location of the charge distribution at a, and in 
the fields atr>b. 

• "Furry's" perspective (a ^ 0, (3 = 0) : 

pa— a p pb poo 111 

The electromagnetic momentum is at the location of the current distribution at b, and in 
the fields at r < a. 

In these symbolic equations the first and fourth integral signs correspond to the integrals in 
the regions (0 ^ r ^ a) and (6 ^ r ^ 00), respectively; and the two integral signs in between 
them correspond to the contributions of the integral over the region (a ^ r ^ 6), which are 
non-zero only when the singularities at the lower and/or upper bounds are included. 

Mathematically, these four possibilities have their origin in the properties of the 
singularities of the fields at a and b. The contributions arising from these singularities are such 
that the integral P = j duj drr'^E x B/An has a unique value which is the same whether 
any singularity is included or not (just like in complex analysis contributions of poles can be 
included or excluded by an appropriate choice of the integration path). 

In physical language, this corresponds to the well-known "wave-particle" duality of field- 
theory, which means that there is as much information in the "fields" than in the "particles" [[8]|, 
a property that in mathematics is systematically exploited in the theory of complex and hyper- 
complex functions. Thus, our result is a beautiful confirmation of Lanczos's intuition, namely 
that electrodynamics can be interpreted as a four-dimensional generalisation of Cauchy's 
theory such that singularities correspond to electrons [14], and is therefore providing another 
example of the remarkable diversity of mathematical and physical pictures which can be 
invoked to interpret electromagnetic phenomena. 

Finally, the reason why the symbolic equations (17.31) and (|7.4I) have been labeled as 
"Faraday's" and "Furry's" is that they suggest an interpretation of the Faraday and Furry 
formulae with regards to the location of the electromagnetic momentum. Indeed, in order to 
integrate over the charge and current distributions as in the these formulae, it is necessary that 
the singularities in the field strengths have been included (either explicitly as with the method 
of reference [3], or implicitly as is usually the case), because otherwise the charge and current 
distributions would be zero. 
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8. Appendix: Poynting's method for finite dipolar system 

In this appendix we calculate the electromagnetic momentum with Poynting's method for the 
double-shell system of section[6l leaving the radial integrals in parametric form with the upper 
or lower integration bounds written as a ± a and b± P, where a <ti a and -C 6 are positive 
infinitesimals. Thus we calculate the momenta separately in the three regions ^ r ^ a, 
a ^ r ^b, and 6 ^ r ^ oo according to equation (11.11) . but leave the radial integrals as given 
by (16:241) . i.e., 

pa— a fb—fS POO 

/ dr {...), / rfr (...), and / dr {...), (8.1) 

Jo Ja+a Jb+13 

This enables us to find out where the electromagnetic momentum is located in the system, as 
well as the respective amounts of the various regular and singular contributions to the total 
momentum. 

As in sections m and m we explain every step in detail. This gives an opportunity to show 
how expressions containing both T and 5 functions can be integrated using simple techniques. 

8.1. Internal region ^ r ^ a 

In the internal region E is the constant electric field (16.81) and B the constant magnetic field 
(16.1 II) . Thus, we have to calculate the vector product E x B where 

Edm{r, r ^ a) = EoT{a - r) - u{Eo ■ u) r5{a - r), (8.2) 
Mr,r^b) =BoT{b-r)-^(^Bo-u{Bo-u))r6{b-r), (8.3) 

where we have used the identity u x (Bq x u) = Bq — u{Bq ■ u) to replace the cross-products 
by a dot-product. In the general case, this means that we have to integrate 6 — 1 = 5 terms 
taking into account that u x u = 0. However, if we avoid the degenerate case a = b and 
specialise to the case 6 > a, we remark that 

T(6-r) = l and S{b - r) = 0, (8.4) 
so that we need to integrate only two terms, namely 

Ex B = EqX BoT{a - r) - u{Eq ■ u) x Bq r6{a - r). (8.5) 
Integrating over angles, this is 

jj dco E X B = EqX Bo(T{a-r) -^r5{a-r)y (8.6) 
-Ppoynting Is thus givcn by the radial integral 

Ppoynting = Eo X Bq dv (^r^T{a - v) - K^6{a - r)) . (8.7) 
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To integrate we notice that 

(-r^T[a — r)\ =r'^T{a — r) r^6{a — r), 

V 3 / 3 



so that the radial integral becomes 

-Ppoynting = EoxBo 

which gives 



dr { -r^T(a — r] 
\ 3 



^'poynting = Eq X Bo^r'T (a - r) 



1 

'3' 







Thus, using (16.131) . i.e., Eq = —p/o.^, we obtain finally 

1 

3 a 



-Ppoynting(0 ^ r ^ o - tt) = ^7^T(a - r] 



Bo X p. 
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(8.8) 



(8.9) 



(8.10) 



(8.11) 



8.2. Intermediate region a ^ r ^ b 



In the intermediate region E is the dipolar electric field (16.21) and B the constant magnetic 
field (16.111) . Thus, we have to calculate the vector product E x B where 



Edmir.r > a) 



P 



T(r — a) 



^\ ' 5{r-a), (8.12) 



Bd{r,r^h) =BoT{b-r)--(^Bo-u{Bo-u))r6{b-r), (8.13) 

where we expressed B^ as in (18.31 ). In the general case, this means that we have to integrate 
9 — 2 = 7 terms taking into account that ux u = 0. However, if we avoid the degenerate case 
a = bvje remark that 

6{b-r)S{r-a) = 0, (8.14) 

so that we need to integrate only six terms, namely 

1 



Ex B 



3u{p- u) X BoT{r — a)T{b — r) 
3 

— ■ u) X BoT(r — a)r6{b — r) 

— px BQT{r — a)T{b — r) 



+ -px [Bq — u{Bq ■ u))T{r — a)rS{b — r) 



u{p ■ u) X Bor6{r — a)T(b — r) 



(8.15) 

(8.16) 
(8.17) 
(8.18) 

(8.19) 



Integrating over angles, this is 



2 y ^ T(r — a)5{b — r) 



-(- 



+ -6{r - a)T{b -r))px Bq. 



(8.20) 



As T(r — a)6{b — r)= 5{b — r) and 5{r — a)T(6 — r) = 5{r — a) when a ^ b,v/e can remove 
the T functions. Thus, -Ppoynting is given by the radial integral 

-/ dr{-6{b-r) + -6{r-a))px Bo, (8.21) 

J a+a 



Poynting 



6 
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which is simply 

^Poyntmg(a + a ^ T ^ h - (3) = - o) - {h - r] 



b-f3 ^ 

Bo X p. (8.22) 

a+a 



8.3. External region 6 ^ r ^ oo 

In the external region E is the dipolar electric field (16.21) and B the dipolar magnetic field 



(16.51) . Thus, we have to calculate the vector product E x B where 

Edm[r,r ^ a) = (^3 — ^ ^jT{r-a} ^ — d{r - a), (8.23) 

Bd{r,r^h) =3 ^3^ - T(r - &) + \ >-5{r-h), (8.24) 

where we have used the identity m x (m x m) = m — u{m ■ u). In the general case, this means 
that we have to integrate 12 — 4 = 8 terms taking into account that u x u = 0. However, if 
we avoid the degenerate case a = b and specialise to 6 > a, we remark that 

T(r - a) = 1 and 6{r - a) = 0, (8.25) 

so we need to integrate only six terms, namely 

E X B = — (^—3u{p ■ u) X fh — 3p X u[m ■ u) + px T(r — b) (8.26) 

H — g- (^+3u X rh{p- u)— pxrh + px u{fh ■ u^6{r — b). (8.27) 
Integrating over angles, this is 

— !! dtoE X B = px m(-^T(r - b) + -^S(r - b)] . (8.28) 
47r JJ \ r° 3r^ / 

-Ppoynting is thus given by the radial integral 

Ppoynting =px dr {v - b) + ^6{r - b)'^ , (8.29) 
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which immediately integrates to 

-Ppoynting = p X m^T{r -b) 
Using (|6.14l) . i.e., fh = 6^i?o/2, we obtain finally 



1 63 

-Ppoyntmg(6 + /5 ^ r ^ Oo) = g^^l'^ ^ ^) 



(8.30) 



Bo X p. (8.31) 

b+l3 
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Errata to reference 

Four lines below equation (2.6) the erroneous expression |r| = sgn(r) should be replaced 
by |r| = r sgn(r); in equation (3.4) the factors dOdO should be replaced by dOdcp; and the 
redundant minus signs in equations (4.4), (5.4), (5.5), (6.7) and (6.8) should be removed, so 
that the charge and current densities read as follows 

47rp„(rO = V ■ i„ = e^6{r) (4.4) 
47rJ,(f) = V X H,{r) = 3^6{r) (5.5) 
47rp,„(r) = V • #d„(f) = 3^^5(r) (6.8) 



